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Abstract 

Energy theory for incompressible Newtonian fluids is, in many cases, capable of producing 
strong absolute stability criteria for steady flows. In those fluids the kinetic energy naturally 
defines a norm in which perturbations decay monotonically in time at sufficiently low (but non- 
zero) Reynolds numbers. There are, however, at least two obstructions to the generalization 
of such methods to Oldroyd-B fluids. One previously recognized problem is the fact that the 
natural energy does not correspond to a proper functional norm on perturbations. Another 
problem, original to this work, is the fact that fluctuations in Oldroyd-B fluids may be subject 
to non-normal amplification at arbitrarily low Reynolds numbers (albeit at sufficiently large 
Weissenberg numbers). Such transient growth, occuring even when the base flow is linearly 
stable, precludes the uniform monotonic decay of any reasonable measure of the disturbance's 
amplitude. 
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1 Introduction 

Linear stability theory can predict sufficient conditions for instability, but some sort of nonlinear 
analysis is necessary to determine stability against finite amplitude perturbations. In many fluid 
dynamic systems the notion of nonlinear energy stability is useful for determining conditions for 
absolute stability ^HJ. Energy stability theory exploits the observation that physically natural 
energy functionals may generate norms that can be used to measure the amplitude of disturbances 
about a base flow. Energy analysis typically leads to sufficient conditions for the norms of arbitrary 
amplitude disturbances to decay monotonically in time, a very strong notion of absolute stability. 
For example, steady states of incompressible Newtonian fluids are very generally energy stable 
at sufficiently low Reynolds numbers. In some cases, notably simple Rayleigh-Benard convection 
problems in the Boussinesq approximation, the critical Rayleigh numbers for energy stability and 
linear instability coincide to establish the supercritical nature of the bifurcation point Moreover, 
energy methods can be generalized to derive rigorous bounds on flow quantities for high Reynolds 
(or Rayleigh) number steady, unsteady and even turbulent flows 0|l3Ej- 
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Elastic and kinetic energy budgets have been used frequently in the past for the linear stability 
analysis of dilute polymer solutions in simple flow geometries. Purely elastic instabilities at small 
Reynolds numbers were analyzed by Joo and Shaqfeh [7j (a review can be found in [S]) for an 
Oldroyd-B fluid in a Taylor-Couette geometry. In the total disturbance energy balance, coupling 
between perturbations of the velocity field and gradients of the polymer stress showed the same 
temporal growth as the kinetic energy of flow and were thus identified as the ones causing the 
instability. Similar energy budget considerations were applied and extended by Sureshkumar for 
linear stability analysis of plane Poiseuille flows at larger Reynolds number [SJ where the impor- 
tance of this coupling was confirmed. For growing Weissenberg numbers the coupling of velocity 
perturbations to gradients of normal viscoelastic stresses was found to dominate the destabilization. 

In this paper we investigate the notion of nonlinear energy stability for Oldroyd-B fluids. To 
set the stage and fix notation, first recall the case of plane Couette flow of a Newtonian fluid. The 
fluid is confined between parallel horizontal plates at y = and y = £, with the lower plate at 
rest and the upper plate moving in the x-direction (i) with velocity U. For simplicity, we consider 
periodic conditions in the streamwise (i) and cross-stream (k) directions. Measuring length in units 
of I and time in units of l/U, the velocity vector field u(x, t) = in + jv + kui of an incompressible 
Newtonian fluid satisfies the Navier-Stokes equations 

-^ + u-Vu + Vp= — Au, Vu = 0, (1) 

at Re 

where Re = Ul/v is the Reynolds number defined in terms of the kinematic viscosity of the fluid 
v. 

The kinetic energy evolution equation is 

d 1 P f 1 f 1 du 

— — I \u(x, t)\ 2 dxdydz = — / — — |Vu| dxdydz + / — — -r— dxdz. (2) 
at 2 J J Re J Reoy y=i 

The negative definite term on the right hand side is the bulk viscous energy dissipation rate, and 
the indefinite term is the power expended by the agent enforcing the boundary condition. 

Plane Couette flow, an exact steady solution in this geometry at every Reynolds number, is 

^Couette = iy- (3) 

If we perturb plane Couette flow at time t = with an initial condition 

u(x,0)=ii/ + (5u(x,0), (4) 
then the energy in the subsequent disturbance <5u(x, t) evolves according to 

— - / \Sulx,t)\ 2 dxdydz = — / (— — I V5u| 2 + 5u5v )dxdydz. (5) 
at 2 j J Re 

Energy stability theory ensures us that there is a positive critical value ReES ~ 83 such that 

Re < ReES guarantees that there exists a number a > such that 

/(— — I V<5u| 2 + 5u5v)dxdydz > — / \5w\ 2 dxdydz . (6) 
Re 2 J 



This in turn implies 



J \5u(x,t)\ 2 dxdydz < J \5u(x,0)\ 2 dxdydz x e~ at (7) 
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so the norm of the perturbation decays monotonically (and exponentially) as t — > oo. 

This is the sense in which plane Couette flow is absolutely stable to arbitrary (finite energy) 
disturbances when Re < Rbes- Plane Couette flow is linearly stable at all Re but finite 
amplitude disturbances may not decay sufficiently high above the energy stability limit. Non-normal 
amplification of small disturbances, a transient growth phenomenon available even in linearly stable 
systems, is recognized as an effective mechanism for stimulating finite-amplitude instabilities 10, 

□2]. 

In the following sections we address the question of whether or not this notion of energy sta- 
bility may be fruitfully generalized to Oldroyd-B fluids. There have been previous analyses in this 
direction for other non-Newtonian constitutive relations |13) . and a priori energy estimates for 
Oldroyd-B fluids 14 , but to our knowledge there has been no definitive determination of whether 
an energy stability principle exists for Oldroyd-B fluids. Here we show by example that no such 
principle is possible in general for Oldroyd-B fluids. By examining particular exact solutions to 
the full nonlinear equations of motion we determine that some perturbations will grow, absorbing 
energy from the base flow, even at arbitrarily small (but nonvanishing) Reynolds numbers. 

2 Energy evolution for Oldroyd-B fluids in the Couette geometry 

In nondimensional form, the equations of motion for an Oldroyd-B fluid in the plane Couette 
geometry 1 are 

-r^ + u- Vu + Vp= — Au + V-r, V-u = 0, (8) 
at Re 

where Re = Ut/v is the Reynolds number. The polymer stress tensor r is 

I) (9) 



ReWi 



where the Weissenberg number Wi = XJXjl is the product of the polymer relaxation time A and 
the imposed rate of strain U/£, the parameter s is the coupling constant proportional to the 
concentration of the polymers in the flow, and the polymer configuration tensor c(x, t) evolves 
according to 

dc 1 

— +u-Vc = c- Vu+(Vu) T -c + — (I-c). (10) 
at Wi 

The component Cj,- of the configuration tensor c is the end-to-end moment (r^j) of the extension of 
the polymers, modeled as linear springs. In this scaling, the equilibrium configuration corresponds 
to an isotropic distribution with unit end-to-end displacements. The elastic energy (density) stored 
in the polymers is proportional to \trc ~ \tr(r) + constant. 

The kinetic energy evolution equation is 

— — l |u(x, t)\ 2 dxdydz = — I — — |Vu| 2 dxdydz + / (tt~^ — YT\2)dxdz— \ r:\7udxdydz. (11) 
at 2 J J Re J y=1 Re Oy J 

The first (negative definite) term on the right hand side is the bulk viscous energy dissipation rate; 
the second (indefinite) term is the power expended by the agent enforcing the boundary condition; 

lr That is, for fluid confined between parallel horizontal plates at y = and y = 1 with the lower plate at rest and 
the upper plate moving in the i-direction (i) with velocity 1 and, for simplicity, periodic conditions in the streamwise 
(i) and cross-stream (k) directions. 
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the third (indefinite) term corresponds to the rate of work done by the polymers on the flow, an 
energy exchange term. The elastic energy stored in the polymers evolves according to 

— - J tr(r)dxdydz = J t : Vudxdydz — J tr(r)dxdydz. (12) 

The two terms on the right hand side are, respectively, the rate of work done by the flow on the 
polymers and the dissipation rate due to polymer relaxation processes. Hence the total energy is 
identified as the sum of the fluid's kinetic energy and the polymers' elastic potential energy, 

£ (t) = J ^[\u(x,t)\ 2 + tr(T)]dxdydz, (13) 

and the total energy evolution is 

- = -J -\V^dxdydz - — J tr(r)dxdydz + J^i.— + ru)dxdz, (14) 

expressing the net balance between dissipation (the first two negative terms on the right hand side) 
and the work done on the fluid by the agency imposing the boundary condition (the indefinite 
surface integral of the stresses) . 

Two distinct obstructions stand in the way of attempts to generalize the property of energy 
stability to Oldroyd-B fluids. The first problem is evident already in (|13|) : while the kinetic 
energy defines a natural norm on the velocity field, the polymer energy does not correspond to 
a norm (or even a metric) on the polymer stress field |14j . Indeed, the physically relevant space 
for the polymer configuration tensor c is the set of positive symmetric matrices, but this does not 
constitute a linear vector space. So strictly speaking the concept of a norm is not meaningful. More 
importantly, though, while it is true that for positive symmetric matrices tr{c) = implies c = 0, 
the trace does not satisfy a triangle inequality necessary for a metric. It is clear, then, that the 
most straightforward generalization of energy stability will not work for Oldroyd-B fluids. 

Nevertheless, various norms or metrics can be defined on perturbations (<5u, 5t) of the variables. 
For example, an I? norm could be defined by 

||<5u, <5t|| 2 = J [5ui5ui + 5Tij5rij]dxdydz, (15) 

Or a functional like 

M(5u,5t) = J [5ui5ui + sj '5Tij5rij]dxdydz (16) 

could be used as an effective measure of the amplitude of perturbations. In any case there is a second 
fundamental — in this case dynamical — obstruction to a principle of energy stability in Oldroyd-B 
fluids: As will be shown in the next section, the problem is that it is possible for perturbations to 
absorb energy from a base flow and, by any reasonable measure, grow at arbitrarily low Re. That 
is, it is not possible to assert in any generality that there is an analog of energy stability in the sense 
that there is a positive critical value ReES s ° that Re < ReES implies that all perturbations decay 
monotonically. This "no-go" result is established by producing specific examples where disturbances 
do not decay monotonically at arbitrarily low Re. 
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3 Some exact solutions in the Couette geometry 



Plane Couette flow ucouette = iy is an exact solution of the Oldroyd-B equations (|8ll0j) in the plane 
shear geometry when augmented with the constant (in space and time) polymer stress tensor 




CCouette = \ Wi 10. (17) 



This flow is linearly stable for upper convected Maxwell fluids ^H], i.e., when s = 0{Re) — > oo. 
Consider perturbations with 5u(x, t) = and 

^cn(i) 5c 12 (t) 0\ 
5c(t) = [ 5c 21 (t) tfpa (t) (18) 
0/ 

where, preserving the symmetry, 5ci 2 = 5c 2 i. That is, we look for solutions where the polymer 
stress tensor is a function of time, but not space, so there is no influence of the polymers back 
onto the flow field. For this class of solutions there are no perturbations of the velocity field or the 
Newtonian stress. 

Inserting this ansatz into the full nonlinear equations of motion produces an exact closed set of 
linear ordinary differential equations for the components of 5c: 

' /<V " L 5 Cll + 25c 12 , (19) 



dt Wi 
d5c\2 1 
dt " ~Wi 

d5c22 1 



5c 12 + 5c 2 2, (20) 
5c 22 . (21) 



dt Wi 
The general solutions, uniform in Re, are 

5cn(i) = [5cn(0) + 25c 12 (0)t + 5c22W 1 ]e- t ' w \ (22) 
Sc 12 (t) = [5c 12 (0) + 5022(0)^-*/™, (23) 
6c 22 (t) = 5c2 2 (0)e' t/Wl . (24) 

These solutions display non-normal transient growth for a variety of initial data as long as Wi 
is not too small. For example consider the zero energy (i.e, trace 5r(0) = 0) perturbation 



5c(0) = a (25) 




where < a < 1. The time-dependent perturbation of the conformation tensor is then 

<t 2 -\ t 0\ 

5c(t) = ae~ t/Wl [ t 10. (26) 
0/ 
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While the perturbations eventually decay exponentially (with algebraic modification), the transient 
dynamics can by highly non-monotonic. For example the polymer-induced stress ~ Ci2(i) overshoots 
its steady value before relaxing. The energy perturbation in the polymers for this solution is 

5£{t) _ tr5c(t) _ at 2 e-*/ m 



Ssteady-state 3 + 2Wi 2 3 + 2Wi 2 ' 

This disturbance initially grows nearly quadratically in time to a maximum at t = 2Wi where 

S£ , 4a Wi 2 

Ssteady-state ' t=2Wi ~ ^3 + 2Wi 2 ' 



(28) 



constituting an Oil) variation in polymer energy before decay sets in. More significantly for our 
purposes is the evolution of a norm-like measure of the perturbation such as the amplitude 



m ysc^jt) I | * (29) 



y/5c tJ (0)5c i3 (0) V 2 

Although A{t) initially decreases, it displays a period of increase when Wi > ^ « 1.043. This 
behavior is plotted for several values of Wi in Figure 1. 

We have found non-normal features in the solutions for any Wi > 1. Consider, for example, 
initial perturbations with 

/0 a 0\ 

<5c(0) = a (30) 
\0 0/ 



with \a\ < \/T + 2Wi 2 — Wi to maintain the positive definiteness of c. The amplitude of this 
disturbance evolves according to 

= ^plpl = e-Wy/TTW, (31) 

and it is easy to confirm that this amplitude exhibits transient amplification so long as Wi > 1. 

A particular example that gives some insight into a physical mechanism for the transient growth 
behavior is the zero energy initial perturbation 

/ -a aWi 0\ 
<Sc(0) = aWi a (32) 
\ 0/ 

where \a\ « 1 corresponds to a (small) uniform rotation of the polymers from their stationary 
configurations. As illustrated in Figure 2, the initially rotated ellipsoid representing the mean poly- 
mer extension is stretched by the background shear well beyond the initial value — by an 0(Wi 2 ) 
factor at large Wi — before reassuming the stationary configuration. During the period of tran- 
sient growth, the polymers absorb energy from the background shear flow faster than the elastic 
relaxation can dissipate it. 
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4 Summary and discussion 



The conventional notion of energy stability of a stationary flow implies monotonic decay of some 
meaningful measure of the magnitude of arbitrary amplitude disturbances. The exact solutions 
developed in the preceeding section show that it is generally not possible to establish monotonic 
decay of arbitrary perturbations of steady flows of Oldroyd-B fluids at finite Reynolds numbers. 
These examples require a nonvanishing (and not arbitrarily small) Weissenberg number so there 
must be some nonvanishing shear in the base flow, but the observed transient amplification of 
perturbations, up to OiWi 2 ) in magnitude, are uniform in the fluid's kinematic viscosity and Re. 
While a base flow may be linearly stable, the 'extra' degrees of freedom in the polymers may allow 
for exact solutions displaying non-normal transient growth. 

We remark that the problems extending energy stability techniques pointed out here are not 
mitigated by using another physically motivated energy functional where a thermodynamic free 
energy takes the place of the mechanical elastic energy alone |16| . That is, a free energy density of 
the form 

F(c) ~ tr{c) - log detc (33) 

includes entropic effects so that the free energy minimum is precisely the equilibrium polymer 
configuration with Cij = 5ij. However this still does not induce a proper metric or norm on 
configurations and because it does not affect the equations of motion, the non-normal growth 
phenomena and the obstructions it presents are still present. 

Finally, after this work was completed we became aware of the analysis of Morozov and van 
Saarloos on an elastic finite-amplitude instability for sufficiently high Weissenberg number 
and arbitrarily low Reynolds number. That work goes beyond the observations presented here by 
explicitly displaying an additional stationary solution to the fully nonlinear 3-d equations. 
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Figure 1: Transient growth of the perturbation amplitude A(t) as a function of time and Weis- 
senberg number as given by Eq. (|29|). 
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(a) steady state configuration 



(b) perturbed configuration at t = 



(c) transient stretch and growth 



(d) relaxation back to steady state 




Figure 2: Sketch of the sequence of deformations of a polymer chain during the non- normal ampli- 
fication phase in a shear flow. 
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